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Abstract. The subject of this article is the Kac equation without cutoff. We 
first show that in the asymptotic of grazing collisions, the Kac equation can 
be approximated by a Fokker-Planck equation. The convergence is uniform in 
time and we give an explicit rate of convergence. Next, we replace the small 
collisions by a small diffusion term in order to approximate the solution of 
the Kac equation and study the resulting error. We finally build a system 
of stochastic particles undergoing collisions and diffusion, that we can easily 
simulate, which approximates the solution of the Kac equation without cutoff. 
We give some estimates on the rate of convergence. 



Mathematics Subject Classification (2000): 82C40, 60K35. 

Keywords: Kinetic Theory, Kac equation, grazing collisions, Particle systems. 

Nicolas Fournier: Laboratoire d'Analysc ct dc Mathematiques Appliquees, CNRS 
UMR 8050, Universite Paris-Est, 61 avenue du General de Gaulle, 94010 Creteil 
Cedex, France 

E-mail address: nicolas.fournier@u-pec.fr 

David Godinho: Laboratoire d'Analyse et de Mathematiques Appliquees, CNRS 
UMR 8050, Universite Paris-Est, 61 avenue du General de Gaulle, 94010 Creteil 
Cedex, France 

E-mail address: david.godinhopereira@gmail.com 

Acknowledgments: The first author was supported during this work by the 
grant from the Agence Nationale de la Recherche with reference ANR-08-BLAN- 



1.1. The model. The spatially homogeneous Boltzmann equation (see Cercignani 
2J, Villani [24]) describes the density ft(v) of particles in a gas, which move with 
velocity v € ]R 3 at time t > 0. The Kac equation is a one-dimensional caricature of 
the Boltzmann equation. It writes 



where t > 0, v S K and where the post-collisional velocities are given by 



0220-01. 



1. Introduction 



(1.1) 




(1.2) 



v' = v cos 9 — sin 9, v' if =v sin 9 + cos 
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The function f} : [— 7r,7r] — {0} — » R+ is an even function called cross section. Each 
pair of particles with velocities v and v* collides to give particles with velocities 
v' and v'^ with a rate proportional to (3(9). See Kac [T3] and Desvillettes [5] for 
more precisions. If we have J Q (3(9)d9 = oo, then there is an infinite number of 
collisions for each particle during any time interval. The case where we assume 
J Q (3(9)d9 < oo (case with cutoff) has been much studied. We will focus here on 
the real physical situation where we only assume f£ 9 2 j3(9)d9 < oo (case without 
cutoff). By analogy with the 3d-Boltzmann equation, we will include the case 
where, for some < v < 2, 

(1.3) (3(9) %° 



We will use in this article Wasserstein distances. Let us recall that for p > 1, if 
/ and g are two probability measures on R with a moment of order p, 

W p (f,g) = inf {E(|C/ - V\ p ) 1/p , U ~ f,V ~ g} , 

where the mfimum is taken over all random variables U with law / and V with 
law g. See e.g. Villani [251 for many details on the subject. In particular, it is 
known that the mfimum is reached : one can build U ~ / and V ~ g such that 
W*(f,g) = E(\U-V\*). 

1.2. Asymptotic of grazing collisions. Assume that there are more and more 
collisions, but that these collisions generate smaller and smaller deviations. For 
example, consider j3 e (6) = jsPy ^rjl|0|< e . Then, we have J* 9 2 f3 e (9)d9 = const 

and J* 9 4 (3 e (9)d9 — > 0. It is known that in this case, the solutions of Boltzmann's 
equation converge to the solution of the Fokker-Planck-Landau equation. To be 
more precise, Degond and Lucquin-Desreux [3] and Desvillettes [3] have shown the 
convergence of the operators (not of the solutions) and Villani [53] has shown some 
compactness results and the convergence of subsequences. The uniqueness results 
of [5] show the true convergence (under some more restrictive assumptions). 

If we denote by (sf)i>o the solution of equation ( |1.1[ ) with cross section /3 e and 
initial condition g^v) — go(v) and if we assume J R v^g^^dv < oo we will show 
that sup te[0 oo) W 2 (gl,gt) < Ce, with (g t )t>o starting from 50 and solving 

(1.4) ^^ = f^W + ^^^' 

where £ := J R v 2 go(v)dv. This limit equation is nonlinear, but the nonlincarity 
appears only through J R v 2 g t (v)dv, which is constant in time. 

A similar result has already been proved by Toscani [5U] with a stronger distance 
but the rate of convergence is not very explicit. We believe that the present rate of 
convergence is optimal. 

1.3. Replacing grazing collisions by a small diffusion term. We come back 



to the Kac equation ( 1.1 ) with fixed cross section (3. Numerically, we must truncate 
small collisions, since they are in infinite number. There are two possibilities. 
One may truncate roughly small collisions by replacing j3 by (3 t (9) = (3(9)l\g\ >e . 



We denote by (/ t c )t>o the solution of (1.1) with this (3, 
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One may replace small collisions by a small diffusion term in the spirit of grazing 
collisions. We denote by (/ t e )oo the solution to 



(1.5) 



where 



(1.6) 



d_ 



ft(v) =b t -0-\vft(v) 



do 



|e|>e 



b f = 



(1- cos 6)j3(9)d6 and £ 



|0|<6 



We will show that sup te[0jT] W 2 (f t ,f£) < Ce(l + \/T) if J R v i f (dv) < oo. Ob- 
serve that when neglecting roughly grazing collisions, we get sup t6 r 0T ] W2(/ t e , ft) < 



Ore 1 -"/ 2 (see Desvillettes-Graham-Meleard [7]) if /3 is as in (TL3J) . We can yet no- 



tice that there is no dependance on v in our result. This is due to the fact that 
the more v is close to 2, the more we neglect small collisions, but the more small 
collisions are well-approximated by the diffusion term. The proof is inspired by [5] . 



1.4. A finite system of stochastic particles. Let j3 be a given cross section and 
fo an initial datum with f R v 4 fo(dv) < oo. We consider a solution (ft)t>o °f (1-1 1. 



For e > fixed, we are going to build a system of n stochastic particles that we 
can simulate with a cost of order Tn f,g, >e (3(8)d8 on [0,T], which is at worst of 

order Te~ 2 n. If we denote by ^™' e the empirical measure associated to this system 
of particles and by /x" the empirical measure associated with a system of n i.i.d. 
particles with law f t , we will show that 



sup E 

[0,T] 



Wi(ft,tf' e )\ < C(l + T) 3 (e 2 + supE[wJ(/ t ,#)l). 

J V [ 0T] L J/ 



Our system of particles is thus as efficient as the system of particles with true i.i.d. 
particles with law ft which are not simulable because of the nonlinearity. If we 
assume that fo has infinitely many moments, we will get 

1 



sup E 

[0,T] 



iyK/ t , M n]<c(i + r)3( e2 + _). 



This system of particles uses the ideas of the previous section : we replace small 
collisions by a small diffusion term, which gives an error of order e. 

In Desvillettes-Graham-Meleard [7], they just cutoff small collisions and they 



E 



get, roughly, something like supj T ] 
A e = f m>e 0(9)d6 w e~ u if /3 is as in (|1.3) 



Wl(futi' e ) 



< G 



T I c 



with 



If we compare this result with our 



result, we can observe the following. 

• In the first term, we get an error of order e 2 instead of e 2 ~ u . It is due to 
the fact that we replace small collisions by a small diffusion term. 

• In the second term, we get a bound which does not depend on e. It is be- 
cause we use a Wasserstein distance which is well-adapted for this study. In 
Desvillettes-Graham-Meleard [7], they give the final result with a Wasser- 
stein distance, but to get this result they use a variation distance. 

• The cost of simulation for the two systems of particles is similar. 
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See also Peyre |16j who gives large deviations estimates for the Boltzmann equa- 
tion for Maxwell molecules and Mischler-Mouhot [TS] who give results of chaos prop- 
agation with quantitative estimates for the Boltzmann equation for hard spheres 
and for Maxwell molecules. 

1.5. Comments. We managed to obtain some bounds uniform in time for the 
asymptotic of grazing collisions. For our two other main results, we tried to limit 
the time dependance. We thus avoid getting bounds with exponential terms. 

The bound we get for E[W r |(/ t , //")] is not very satisfactory. A priori, it is of 
orde r n^ 1 / 2 ) - (if the initial condition has infinitely many moments, see Lemma 
A. 4) which gives a bound for E[W 2 (ft, £t?)] °f order n - ^ 1 / 4 ' - . We expected to get 



a bound of order n" 1 / 2 as in the central limit theorem, but we cannot get it. See 
Peyre [T5] for example to get more details. It seems to be the only defect of W 2 for 
this study. 



Assuming that J* 0/3(0) d0 < oo (e.g. if we assume (1.3 1 with v € (0,1)), we 
get a bound for E[Wi (ft, which is of order e + n -1 / 2 but with an exponential 

dependance in time. If J* 1 /3(0)dO < oo (e.g. if v < 7) for some 7 G (1,2), we 
also study E\W](f t , £>')]. 

In a future work, we will apply the same kind of methods to the homogeneous 
Boltzmann equation. We hope to get some results which will probably be much 
less optimal. 

Our proofs use probabilistic methods, which was initiated in the famous paper 
of Tanaka [T5], and used in Desvillettes-Graham-Meleard [7]. We will also use a 
result of Rio [TB] , which gives some very precise rate of convergence for the standard 
central limit theorem in Wasserstein distance. 

1.6. Plan of the paper. In the next section, we will state more precisely our 
three main results. In Section 3, we will give a probabilistic interpretation of the 
three equations. Sections 4, 5 and 6 are devoted to the proofs of our main results. 
Some numerical illustrations will be given in Section 7. At the end of the paper, 
we will give an appendix with some results about the Wasserstein distance between 
a compensated Poisson integral and a centered Gaussian law with same variance, 
the rate of convergence of an empirical measure using Wasserstein distances, the 



moments of the solution to ( 1.1 ) and the well-posedness for a certain kind of P.D.E.s. 



2. Results 

2.1. Weak solutions. Let f3 be a cross section satisfying 
(2.1) / 6 2 f3(0)d0 < 00. 



For k > 0, we denote by 'Pfc(K) the set of probability measures on R admitting 
a moment of order k and by C 2 (M.) the space of real bounded functions which 
are in C 2 (1R) with first and second derivatives bounded. We say that a family of 
probability measures (f t )t>o is in Lf^ c ([0, oo),V 2 (R)) if sup [0T] j R v 2 f t (dv) < 00 
for all T.Ifipe C 2 (R) and (v, v,) £ R 2 , we set 



Kg(v, u*) = / tp(vcos0 — v*sm0) — (p(v) — (v(cos0 — I) — v*smO)(p'(v) (3(0)d0 
(2.2) -bvif'(v), 
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with 
(2.3) 



b = 



(1- cos 9)/3(8)d9. 



If J* 6/3(O)d0 < oo, then one easily checks, using that (3 is even, that 



(2.4) 



K?(v,v*) 



ip(v cos 9 — i>* sin 6) — <p(v) j3(9)d9 



We now define precisely the notion of solutions that we will use. 



Definition 2.1. Consider a cross section j3 satisfying (2.1). 



(1) We say that (f t )t>o € i™ c ([0, oo),^^)) sofces i/ for any ip in 

C£(R), any t > 0, 



(2.5) 



ip(v)f t (dv) = / (p(v)f (dv) + 



Kp(v, v*)f s (dv)f s (dv*)ds. 



(2) We say thay (g t )t>o & L uS c (i Q > °°)> ^(K)) sofoes i/ /or any y> in 

C b 2 (M), any f > 0, 



(2.6) 



tp(v)g t (dv) 



where £ := J v 2 go(dv). 



(p(v)g (dv) 



'-£ 



<p (v)g s (dv)ds 



2 Jo 



vip (v)g s (dv)ds, 



(3) For e G (0,1) /wed, we say t/iai (/ t e ) t > G i^ c ([0, oo), P 2 (M)) so/ves |X5) 
if for any ip in C 2 {M), any t > 0, 



(2.7) 

where 



<p(v)ft(dv) 



v(p'(v)fg(dv)ds 



+ £b f 



^"(v)fl{dv)ds 



Kl(v,v*)r s (dv)fl(dv*)ds, 



(2.8) p e (6) = j3(9)± m>e , b e = {1- cos 9) P(9)d9 and £ = / v 2 f^{dv). 

J\8\<e Jr 

Observe that all the terms in the above equations arc wcll-dcfincd. For exam- 
ple in (2.5 1, the last term is well-defined because for <p G Cj(R), \Kf{v, u*)| < 



Proposition 2.2. Let /o, 50 /o oe * n ^(K) and Zei /3 satisfy (2.1). There 
is existence and uniqueness of solutions (ft)t>o> (.9t)t>o an d (ft )t>o to equations 



1.1), {1.4) an d (1.5) starting from fo, go and /q respectively, in the sense of 



Definition \2.1\ Furthermore, we have energy conservation: for any t > 



(2.9) 

and 

(2.10) 



v 2 ft(dv) 



v 2 fo(dv), / v 2 g t (dv) 



v 2 go(dv) 



v 2 ft(dv) 



v 2 r {dv). 
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F or th e proof of the pr eviou s result, one can see Toscani-Villani [21] for 

with a = i£, b = — \ and q = r = 0. For 



For (1.4), use Proposition 



use Proposition A. 6 

by r(t,v,v*,A) 



A.6 



with a = Eb t , b 
j! 7r l|0|>elA(u(cos 9 



1.1 
1.5 



-6 £ , q = and r(t,v,v*,dh) defined 
1) — v* sin 9)(3(9)d9 for all Borel subset 
A of E, which indeed satisfies sup t v ^ r(i, u, u t , E) = J*^ li0i >e /3(0)d# < oo and 
su Pt > J R (/i 2 +2«ft)r(t,«,«„dfc) = L sin 2 9 P(9)d9(v 2 + v 2 ) = C{v 2 +v 2 ). To 



get energy conservation, it suffices to apply (2.5), (2.6) and (2.7) with ip(v) 



2.2. Asymptotic of grazing collisions. Our first main result is the following. 

Theorem 2.3. Let go € ^(E) and let (/3 e ) e g(o.i) be a family of cross sections 
verifying f_J 2 fi t (9)d9 = 1 and J^9 i l3 e (9)d9 ^> 0. For e € (0, 1), let (g e t ) t >o be 



the solution of 
the solution of 



l.ly with go for initial datum and /3 e for cross section. If (gt)t>o is 
XjP with the same go for initial datum, then for all e € (0, 1), 



sup W 2 (gt,g t ) <C 

te[o,oo) 



f M v 4 go(dv) 



9 4 p e (9)d9, 



where C is a universal constant. 



This convergence result was already known (see for example Toscani [20]), but 
we get here an explicit and probably optimal rate of convergence, which, to our 
knowledge, had never been done so far. 

Remark 2.4. If we consider a cross section /3 with J^9 2 f3(9)d9 — 1 and if for 
any ee (0,1), we set ft(0) = ^(f )t\e\ <e , then 9^ t (9)d9 < e 2 . 

2.3. Error when we replace the small collisions by a small diffusion term. 

Let us explain briefly why (2.7) approximates (2.5): consider a cross section j3 



satisfying (2.1) and ip € C£ 



Write, using that (3 is even, 



<p(vcos9 — sin6>) — tp(v) — (v(cos9 — 1) — sin 6) ip'(v) j3{9)d9 



+ K^(v,v*) - b e vtp'(v) 
v(cos 9 — 1) — sin ( 

e 

sin 2 9/3{9)d9 



(3{9)d9 



+ Kt (v, v*) - b e v(p'(v) 



2 



2 



2 ' "A 

(cos9~l) 2 /3(9)d8 



+ K^(v,v^) - b e vip'(v) 
« v"{v)v% +Kl(v,v*) - b e v<p'(v). 

We decided to neglect the second term in the fourth line of this approximate equal- 
ity, since it is much smaller than the other terms, because ft e ( cos @ — l) 2 (3(9)d9 < 

In order to obtain an equa- 



fie 8 4 P(8)dQ < e 2 Il e 2 li{9)d9 a e 2 £ e sir 
tion preserving the kinetic energy, we replaced | f* sin 2 9(3(9) d9 by 6 e (both are 
approximately equal to \ j e _ e 9 2 (3(9)d9). Our second main result is the following. 
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Theorem 2.5. Let fo 6 "P^R) and £ef ft be a cross section satisfying (2.1). For 
e € (0,1), we consider (ft)t>o an d (ft)t>o solutions of (1.1) and (1.5) respectively, 
both starting from fo. Then for any T > 0, any e £ (0, 1), we have 

1 



sup W^ftJl) <CVmin(l + T, 
te[o,T] 



/I 



|e|<6' 



2 /3(0)d0 



where C depends only on J R v 2 fo(dv) , J R v 4 fo{dv) and on J™ 9 2 (3(9)d9. 



We can observe that we are not so far to get a bound uniform in time for e 2 (we 
do not have exponential bounds). 



Remark 2.6. If (3 is as in 



1.3), we get a bound in C min ^e 2 (l + T), e"^ . 



2.4. System of particles. Let fo be a probability measure on K and let (3 be a 
cross section satisfying (2.1 ). We fix an integer n and we consider: 

• a family of i.i.d. random variables {Yo)i&{i,...,n} with law fo, 

• a family of i.i.d. Poisson measures (N l ' n (dsd9dj)) il£ ^ lt n j on [0,oo) x 
[— tt,tt] x {l,...,n} with intensity measure ds)3(6)d0^ X^fc=i ^k(dj), 

• a family of i.i.d. Brownian motions (i?|) t > 0; ...,«}■ 

For e g (0, 1), we consider (V^' n ' e ) t >o, ie{i,...,n} solution of the following system of 

SDEs: fori = 1,..., n, for alU>0,~ ' 

(2.11) 



J\e\>eJje{l,...,n} 

t 



(cos 9- l)V^' e - sin0F/i n,e N l - n (dsd9dj) 



b e f V*' n > e ds + ^/2£b~ e B l t 
Jo 



't 



where £ — J R v 2 f (dv) and b e — Li <e (l — cos9)j3(9)d9. 

The quantity V^' n ' e has to be thought as the velocity of the i-th particle at 
time t. The behavior of (V/'"' £ )t>o is the following: after an exponential time 
r with parameter A e = Jj e i >e (3(9)d9, it collides with another particle labelled j 

chosen at random and then we set Vp n ' e = (cos Q)V^'((( :e — (sin 0)V^?l"' c , where O is 
A~ 1 / 3(0)l|gi> e d6'-distributed. Between two jumps, V l ' n ' e behaves like an Ornstein- 
Uhlenbeck process 

- b e J V*' n - e du + ^2Eb e (B\ - Bl). 
We can solve explicitly this last SDE and we get 

yi,n,e _ yi,n,e e -b4t-s) + y ^£ ~ e -bAt-s) f 

J s 

Hence the strong existence and uniqueness of a solution (V/' n ' e ) t >o, ie{i....n} to 
(2.11) is straightforward. 

We can observe that to simulate our system of particles on [0, T], we need to 
simulate in mean nT J\g\ >t P(9)d9 jumps. We thus have a cost of simulation of 
order nT f^ e ^ >e [3(9)d9. The fact that we can explicitely solve the previous SDE is 
fundamental in order to have such a cost of simulation. 
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Theorem 2.7. Let /q G V4 (R) and Zet f3 be a cross section satisfying (2.1). We 
consider (ft)t>o solution to the Kac equation (1.1) starting from /q. For n G N* 
and e G (0,1), we consider the solution (V t l: "' e ) t >o i to (2.11). We set 

/i™' e = i <5yi,n, e . XVien /or any T > 0, any n > 2 and any e G (0, 1), we /iawe 



sup E 

te[o,r] 



W?(/t, a4" £ )1 < C(l + T) 3 (e 2 + sup E\wi(f t , tf) 

V [0,T] L 



where C depends only on f^v 2 fo(dv) , f K v 4 fo(dv) and on J* 9 2 (3(9)dQ , and where 
for all t > 0, /it" = i 5y», where (V(f)i£{i ,...«.} * s a family of i.i.d. particles with 
law f t . 



Applying Lemma |A.4| of the appendix we will deduce the following consequence: 



Corollary 2.8. Under the same assumptions and notation as in Theorem \2.1[ if 

/o has a moment of order p > 4 with p even, then for all T > 0, all n > 2 and all 
e 6(0,1), 



sup E[^ 2 2 (/ t ,^' e )l <C(1 + T) 3 ( £ 2 + ^: 



te[o,T] 



_2 



where C depends only on p, fo and L 9 2 (3(9)d9 . 

We end this section with a result using another Wasserstein distance. 

Proposition 2.9. Under the same assumptions and notation as in Theorem \2. r l\ 

if the cross section (3 satisfies the stronger assumption 9(3(9)d9 < 00, then for 
all T > 0, alln>2 and all e G (0, 1), 



sup E 

te[o,r] 



where Ct depends only on T, J M v 4 fo(dv) and on 9/3(9) d9. 

We thus have a better dependence in n, but we get exponential bounds in time. 

3. Probabilistic interpretation of the equations 

This section is strongly inspired by Tanaka [19 and Desvillettes-Graham-Meleard 
[TJ. Until the end of the article, (O, J 7 , (.F t ) t >o,P) will designate a Polish filtered 
probability space satisfying the usual conditions. Such a space is Borel isomorphic 
to the Lebesgue space ([0, 1],B([0, l]),da) which we will use as an auxiliary space. 
To be as clear as possible, we will use the notation E for the expectation and C for 
the law of a random variable or process defined on (f2, J 7 , P), and we will use the 
notation E a and C a for the expectation and law of random variables or processes on 
([0, 1], B([0, 1]), da). The processes on ([0, 1],B([0, 1]), da) will be called a-processes. 

We say that a R- valued process (Vt)t>o is a L 2 -process if it is cadlag, adapted 
and if E(sup [0 T1 V 2 ) < 00 for all T > 0. 



Now, we introduce a nonlinear stochastic differential equation linked with ( 1.1 1 



Proposition 3.1. Let ft be a cross section satisfying (2.1\). Let fo G Vz^i) and let 



(ft)t>o be the solution to (1.1) starting from fo. Consider any a-process (Wt)t>o 
such that C a (W t ) = ft for all t > 0. Let also N be a (F t )t>o-Poisson measure on 
[0, 00) x [0, 1] x [— 7r, 7r] with intensity measure dsda(3(9)d9 , and Vq a J~q -measurable 
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random variable with law fo- Then there exists a unique L 2 -process (Vt)t>o such 
that for all t > 0, 

(3.1) 

V t = V - 



o Jo 



(cos0-l)F s _ - sin 6W S _ (a) N(dsdad9) 



b / V s ds, 
Jo 



with b given by (2.3). Furthermore, C(Vt) — ft for all t > 0. 



Proof. While stated in a slightly different way, this result is almost contained 
in Desvillettes-Graham-Meleard [7J Theorem 3.4]. See the proof of Proposition 3.3 
below for similar arguments. □ 



Let us now write down a probabilistic interpretation of (1.4 1. 



Proposition 3.2. Let go G and se t £ = J^ v2 9o(dv). Consider a To- 

measurable random variable Yq with law go and a (jFt)t>o-Brownian motion (B t )t>o- 
Then there exists a unique L 2 '-process (Y t )t>o such that for all t > 0, 



(3.2) 



Y t =Y - 



1 



Y,ds 



f£Bt 



Furthermore, C{Y t ) = g t for allt > 7 where (gt)t>o * s ^e unique solution to (1-4-) 



Proof. The existence and uniqueness of Y is classical since (3.2 1 is a S.D.E. 
with Lipschitz coefficients. By Ito's formula, we have for any ip G Cj(M.) 



<p(Y t ) = <p(Y ) 



<p'(Y s )(--Y s ds 



>£dB s 



£ 



ip"(Y s )ds. 



Taking expectations and setting \x t = £(Yt), we get for any ip G Cf (R) 



<p(v)(i t (dv) = / ip{v)g (dv 



vtp' (v)/j, s (dv)ds + 



ip"(v)n s (dv) 



Thus (/i()t>o solves (1.4 1 in the sense of Definition 2.1 We get (fj, t ) 



by uniqueness (see Proposition 2.2 1 



>t>o 



(gt)t>o 
□ 



It remains to give a probabilistic interpretation of (1.5 1. 



Proposition 3.3. Let e G (0, 1) be fixed. Consider a cross-section ft satisfying 
\2. a probability measure /g G ^(K), and the corresponding unique solution 



(ft)t>o t° (1-5). Consider any a-process (W t e )t>o such that for all t > 0, C a (W^) 



/ t e . Let V e be a !Fo -measurable random variable with law /q 7 let N be a (J r t)t>o- 
Poisson measure on [0, 00) x [0, 1] x [— 7r,7r] with intensity measure dsda/3(9)d9 and 
let (B t )t>o be a (!Ft)t>o-Brownian motion independent of N. Then there exists a 
unique L 2 -process (V t e )t>o such that for all t > 0, 



(3.3) 



J\e\>eJo 
t 



(cos6» - l)V s e _ - sin6W s £ _(a) N(d9dads) 

b e f V s e ds + y/2£b e B t , 
Jo 



with b e defined in (2.8). Furthermore, C(V t e ) — / t e for all t > 0. 
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Proof. See Ikeda-Watanabe [TTJ Theorem 9.1] for existence and uniqueness 
of (V t e ) t >o'- (3.3 1 is a classical jumping S.D.E. with Lipschitz coefficients. Let 
<p G C£(M.). By Ito's formula for jump processes (see e.g. Ikeda-Watanabe [TTJ 
Theorem 5.1]), we have 



<p{Vf) =cp(V c 



o)+ / f'(V:){-b £ V:ds + y/2£b~ e dB s ) + £b e f <p"(V s e )ds 
Jo Jo 



o J\e\>eJo 



ip(cos9V s e _ - sin6W s e _(a)) - (p(V s e _) N{d9dads) 



Taking expectations and setting //| = C(V t e ), we get for any ip € C 2 (R), using that 

ip{v)nl(dv) = [ <p(v)f£(dv) - b e [ [ vLp'(v)^ s {dv)ds 
Jr Jo Jr 

Sb, I I <p"(v)fi e s (dv)ds 

(p(vcos9 — v* sin9) — tp(v) /3(9)d9fi e s (dv)f^(dv^)ds. 



But (ft)t>o solves the same equation since it solves (1.5) in the sense of Definition 
2.1 Since (/|) s >o is given, this equation is linear and we have uniqueness of the 

Sb e , b 



solution. Indeed, we use Proposition |A.6| with a = Eb t , b = — b e , r = and 
q(t, v, A) — J R l\ e \ >e tA(v(cos9 ~ 1) — sin 6*) f%(dv*)P(9)d9 for all Borel subset 
A C R, which satisfies sup t v q(t, v, R) = J, g , f3(9)d9 < oo and sup t>0 / R (/i 2 + 



2vh)q(t,v,dh) = J m>e sm 2 9p(9)d0(f R vlfo(dv« 
get {i4)t>a = (ft )t>o- 



< C(l + v 2 ). Finally, we 
□ 



4. The Grazing collisions limit 
We consider a family of cross sections (/3 e )ee(o,i) with J*^9 2 {3 e (9)d9 = 1 and 
f%0 A p t (9)d9 ^ 0. Let g Q G V^R). For any e € (0,1), we consider (g t £ ) t > the 



unique solution of (1.1) with cross section j3 t starting from go. We also consider 
(<?t)t>o t ne unique solution of (1.4 1 starting from g . For e G (0,1), we consider 



a J-o-measurable random variable Vq with law go, and a (J r t)t>o-Poisson measure 
N e on [0,oo) x [0,1] x [— tt,tt] with intensity measure dsdaf3 e (9)d9. We also con- 
sider an a-process (Wf) t >o such that £ a (W t e ) = <?| for all f > 0. Let (B t )t>o 
be a (J r t)t>o-Brownian motion. We consider (V t e )t>o and {Yt)t>o solutions of the 
following S.D.E.s 



V* = Vo 



o Jo 



(cos 9 - l)V s e _ - sin6W s e _(a) N e (dsdad9) - b e [ V s £ ds, 

Jo 



where b t = f\(\ - cos0)/3 £ (0)d0 and 



(4.1) 



Y t = V -- / Y 8 ds + V£B t . 



Theorem |2.3| is a corollary of the following statement. 
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Theorem 4.1. For any t > and any e S (0, 1), we can couple the Poisson measure 
N t and the Brownian motion B in such a way that 

£f^(l - cos 9) 2 f3 e (9)d9 (E(V 4 ) + 3£ 2 ) sin 4 9 [3 e (9)d9 

T" G 



E[(U t £ - Y t f] <4 



2b, 



+ S(\ ln(26 e )| 2 + + l) \2b t - 1| + 2| 7e - 1|) 

where b e = J^(l - cos 9) fi e (9)d9, 7e = sin 2 9f3 e (9)d9, £ = E[U 2 ] and C is a 
universal constant. 

Let us insist on the fact that the coupling between N e and B depends on t. 
Assuming for a moment that this result holds true, we can prove Theorem |2.3| 



Proof of Theorem 



3.1 



2.3 



and C(Y t ) — g t by Proposition 



First reca lling that C(Vf) = g\ by Proposition 
" we have W 2 (gl g t ) < E[(V t e - Y t ) 2 }. 



3.2 



If f* n 9 i (3 e (6)de > 1, we have E[(V t £ - Y t ) 2 ] < 2E[(V t e ) 2 } + 2E[Y t 2 } = 4£ < 
A£ f* 7r 9 4 j3 e (9)d6. We now suppose that 9 4 (3 e {9)d9 < 1. Using the Taylor- 
Lagrange inequality, we have 

|l-cos0| <9 2 /2, |2(l-cos6»)-6» 2 | < 6> 4 /12 
| sin<9| < \9\ and | sin 2 9 - 8 2 \ < 9 4 /3. 

Using these inequalities, we get 



(1- cos 9) 2 (3 e (9)d9< / (9 4 /4)f3 e (9)d9 



sin 4 9(3 e (9)d9 < 



i l3 e {9)d9, 



and, recalling that 8 2 p e (9)d9 = 1, 



(9 4 / 12) (3 c (8)d9, | 7e - 1| < / (8 4 /3)l3 e (9)d9. 



\2b e -l\ < 



Since J* 7i 9 4 /3 e (9)d9 < 1 by assumption, we have 2b e > 11/12, 2/3 < 7e < 4/3 and 
2b e - 1 £ [-1/12, 1/12] which allows us to write | ln(26 e )| 2 = | ln(l + (2b e - 1))| 2 < 
4|2& £ -1| 2 . WethusgetE[(U t e -F t ) 2 ] < C(£+^p-) J^9 4 /3 e (9)d9, which concludes 
the proof, since £ 2 — E[U 2 ] 2 < E[V 4 ] by the Cauchy-Schwarz inequality. □ 



It remains to prove Theorem 4.1 Let us start with the following lemma. 
Lemma 4.2. For e G (0, 1) 7 let Y t be the unique solution of 



(4.2) 



Yt =V -b e Y^ds 



JO J-tt 



sin 9Wl_ (a)N e (dsdad9) 



(since W e is a given a-process, this is a classical S.D.E. with Lipschitz coefficients). 
Then for all t > 0, 



E 



t - Yt 



e\2 



< 



£ J^(l- cos 8) 2 (3 e (0)d9 



2b t 
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Proof. Observing that 

ft />1 />7T i>t 

V t e -Y f e = / / (cos9-l)V s e _N t {dsdad9)-b e / (V s e -Y s £ )ds 
Jo JO J--K Jo 

we get by Ito's formula 

ft /-l pit r 



JO J -it 



< ^2 



v?_-y;_ + (coB(9-i)K, e _ 

- (V/_ - F/_) 2 1 N £ (dsdad9) 

fo fo I- ^ K ~ ^ + (C ° S ^ " 1)Vs ^ 2 " (V * " 

- 2(V s e - F/)(cos 6> - 1)1*7 dsda(3 e (6)d9 

2b t [\v:-Y:fds 
Jo 



Ml 



t pit pt 

{cos6-l) 2 (V s e ) 2 dsl3 e (9)d6-2b e / (1£ - Y s e ) 2 ds, 

9 J-i JO 



where M t c is a martingale with mean 0. So using that E[(V t e ) 2 ] = J R v 2 g^(dv) 
J R v 2 go(dv) = £ for all t > by (2.9 1, we have 



E[(V< - n) 2 } = [ r (1 - cos9) 2 El(V s <) 2 }p e (e)d6ds - 2b e [ E[(^ e - Y£) 2 ]ds 

JO J -it JO 

= £t[ (1 - cos 9) 2 /3 e (9)d9 - 2b e f E[(V s e -Y s e ) 2 ]ds. 

J -IT JO 

Differentiating this equality with respect to t, we find an O.D.E. that can be solved 
explicitly. This gives 



The conclusion follows. 



£ J^(l -cos 9) 2 (3 e (9)d9 



(l~e- 2M ). 



2b e 



□ 



In the following lemma, using Corollary A. 2 we will find a suitable coupling 
between our Poisson measure iV e and our Brownian motion B. 

Lemma 4.3. Let Y e be the unique solution of 



(4.3) 



V = Vo 



Y^ds + V£ le B t 



We consider the process Y e defined in Lemma 4-2 For any e G (0, 1) and for each 
t > 0, we can couple the Poisson measure and the Brownian motion B in such 
a way that 

i (E(V 4 )+3£ 2 )J^ S hx i 9p e (9)d9 



E[(Y/ - F/) 2 ] < c 



where C is a universal constant and 7 e = sin 2 9(3 e {9)d9. 
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Observe that for each t we need a suitable coupling. We are not able to find a 
coupling working simultaneously for all values of t. 

Proof. Applying Ito's formula, we get Y t e e bet = Vo + y/Sje Jl eheS dB s and 
Y t e e b ' t =V - Jo fo e Ks sin 9W<_ (a) N e (dsdad6). We observe that the random 
variable \f£^u J e b ' s dB s follows a centered normal law with variance £j e L e 2b ' s 
which is equal to /„' £ j\ e 2b ' s sin 2 9(W^(a)) 2 (3 e (9)d6dads because C a {W e s ) = g% 
and due to ( |2.9[ ). So using Corollary A. 2 we get 

Jo Jo /-» e4M sin ' 



^(£(y/e^),£(y t e e ^))<Co 



C 



£ lt Jo e2Md « 
/ 4 e 4 ^E Q ((F^) 4 )d S Z« sm 4 0/3 6 (0)d0 



Using Lemma A.5| since C(Wf) = g e s and since c/ e solves dTTTj l (with the cross section 
P e ), we deduce E Q [(W^ S £ ) 2 ] < J R v 4 g (dv) + 3£ 2 = E[V 4 ] + 3£ 2 for all s > 0. Hence, 

using that 4^ 



Jo e ib ' s d° 
fields 



< e 



2b,- 1 



we have 



W 2 {C{Y?e b ^),C{Y?e b ^))<Cv 
Consequently, 



(E(F 4 ) + 3£ 2 ) sin 4 6(3 t (6)d0t 
ft 



2M 



^|(£(y/),£(y/))<c7o 



(E(U 4 ) + 3£ 2 ) sin 4 0/3 6 (0)d0 



To conclude, it suffices to take N e and B in such a way that E[(Y t e — Y t e ) 2 ] = 
W$(£(Y t % £(¥*)). ' □ 

Let us now give the last lemma needed to prove Theorem |4.1| 



Lemma 4.4. Consider the unique solutions Y and Y e to (4-1) and (4-3) respec- 
tively, driven by the same Brownian motion B. Then for all t > fixed and for all 
£6(0,1), 

Ep7 - Y t ) 2 ] <S(\ ln(26 £ )| 2 + + l) |26 e - 1| + 2| 7e - 1|) . 

Proof. We have Y t c = Vbe~ M + ^/£^e~ brt Jl e b - s dB s and Y t = Vbe"*/ 2 + 

■\f£e~ t l' 2 f Q e s l 2 dB s as in the proof of Lemma 4.3 Since B and Vq are independent, 
we have 



E[(Y t ^Y t ) 2 ]=E(V 2 )(e 



-b.t _ e -t/2\2 



E 

-M 
-6„t 



- t/2 ) 2 + f / 



s / 2 ) 2 ds. 
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We set h(t) = (e bet — e */ 2 ) 2 . The function h reaches its maximum at to = 
ln(26 e )/(6 £ -l/2). Moreover, \h{t )\ = | e -M _ e -V 2 |2 < \ be t -t /2\ 2 = \t \ 2 \b e - 
1/2| 2 = |ln(26 e )| 2 . Next, 



(Vl e e- b < s - e- s ' 2 fds < / {^f t e- b ' s -e- s ' 2 ) 2 d S 



1 



26 £ 6 £ + l/2 

7 £ (6 £ + 1/2) + 2fc £ (fc £ + 1/2) - A^/%b e 



< 



26 £ (6 £ + l/2) 
7e (l/2 - b e ) + 2b £ ( le + b £ + 1/2 - 2^77) 

-(l-26 £ )+2[(,/77-l) 2 + (fo £ -1/2) 

< (^ + l)|26 £ -l|+2| 7 ,-l|, 

the last inequality coming from (s/x — y^y) 2 < \x — y\. The lemma is proved. 



1 

2b, 



□ 



We can now conclude this section. 



Proof of Theorem 4.1. For e G (0, 1) and t > fixed, we take the Poisson 
measure N e and the Brownian motion B as in Lemma 14.31 and we consider the 
processes V e , Y, Y e and Y e build with this N e and this B. Then, writing 



E[(v? - Y t y\ < 4 E[(v t e - r/) 2 ] + E[(y t e - r/)1 + E[(r t £ - y t ) 



and using Lemmas 4.2 4.3 and 4.4 we immediately conclude. 



□ 



5. Cutoff approximation with diffusion 



The whole section is dedicated to the proof of Theorem 2.5 
and let (3 be a cross section satisfying (2.1 ). We fix e £ (0, 1) 
solutions (ft)t>o an d (ft)t>o to (1.1) and (1.5) respectively, both starting from f 



Let thus f E "P 4 (M) 
and we consider the 



We will proceed as follows. We fix some to > for the whole proof. We will 
build some solutions (V t ) t >o and (V t e ) t >o to (3.1) and (3.3), both starting from 
some initial value Vq with law / , coupled in such a way that E[(V to — V t 
small as possible. 



t Q ) 2 ] is as 



We divide the proof into five steps. In the first step, we introduce the (suitably 
coupled) processes (V t )t>o, {Vt)t>o as well as an intermediate process (V t £ ) t >o- 
In Step 2, we upperbound E[(V^ — V^) 2 ]- Step 3 is dedicated to the study of 
E[(y/ o - V U) ) 2 ]. In Step 4, we show that E[(F t e - V t \){V t \ - V to )} = 0. We conclude 
in Step 5. 
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In the whole section, we will use the notation 
(5.1) b e = [ (1 - cos 6) (3(9)de, c € = 2b e + [ sin 2 8f3(9)d9, 



d r = 



[ sin 2 9/3(9)d9 and J e = f (1 - cos 9) 2 p(9)d9. 
J\e\>e J\e\<€ 



l\8\>e J\e\< e 
Step 1: the coupling. 

- Let (Oj, jF l , (Fl)t>o, Pi), i = 1,2, be two Polish filtered probability spaces 
satisfying the usual conditions and consider the following filtered probability space 
(n,T, (J" t )t>o,P) = (Oi x fi 3 , T 1 ® J" 2 , (J-i ® J?) t >o,Pi ® Pa). We denote by E 
the expectation under P and by E, the expectation under Pj. 

- On J 71 , (^^t^OiPi), we consider a /o-distributed random variable V J'q- 
measurable, as well as a (J r t 1 ) t >o-Poisson measure iV|g|> e on [0, oo) x [0, 1] x [— tt, n] 
with intensity measure dsda/3(6)liQi >e d9. We set 

ft pi t* 7T 

(5.2) X t= / / (cos0- l)N m > e (dsdad9) - b e t, 

JO JO J -TT 

We consider the Doleans-Dade exponential of X, see Jacod-Shiryaev jT3l Theorem 
4.61], defined by 

(5.3) Z t = l+ f Z s _dX s . 

Jo 

There holds 

(5.4) Z t = e Xt JJ(1 + AI 8 )e" AXs = e" M [J cos^l^t, 

s<i i>l 

where (Tj, 0j, «j)j>i are the marks of the Poisson measure iV|g|> e . Observe that 

a.s., Z t 7^ Vt > 0, because (3(9)d9 does not give weight to j — |, ~ j. 

Of course, the processes (X t ) t >o and (Zt)t>o depend on e but we do not write 
this dependence in order to lighten notations. 

- For each t > 0, we consider some a-random variables W t and W£ with respective 
laws ft and fl verifying 

(5.5) Wi(f t JZ)=E a [(W t -WZ) 2 }. 

- Recall that t > is fixed. Fix also U3\ G fix. On (ri 2 , J" 2 , (J" 2 ) t >o, P2), we 
consider a (.F 2 )t>o-Poisson measure ^ r |e| <e on [0, 00) x [0, 1] x [— vr, tt] with intensity 
measure dsdafi(9)l\g\ <e d9 and a Brownian motion (B^ 1 ) t >o (we do not write the 
dependence in io and e) such that: 

W 2 2 K 1 5 < 1 )=E 2 [( f° f /^Z to Z7i)( Wl )sin(-0)^_(a)^f <e ( W2 ,d S ^) 



^0 J -It 

2 



(5.6) -jf v^W^to^" 1 )^!)^ 1 ^))' 
where 

(5.7) ^=C 2 (f° ^ y (Z t0 Z7i)(«i)sin(-ff)W a _(a)^ <e (t^,d«dctt»)] 
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(5.8) 



C-2 



( p v / 2^(^ ^- 1 )(^i)^r(a; 2 )). 



Here again we do not write the dependence in e of /i^ 1 and . 

- For (uii,ui2) G we can now set N(u>i, 0J2) — N\ e \ >e (uii) + N^ <<l (uj 2 ) and 
(B t (u}i, W2))t>o = (-St Jl ( w 2))t>o- Clearly, as random objects on (SI, J 7 , T t , P), the 
process (B t )t>o is a (J^t^o-Brownian motion and iV is a (J r t ) t > -Poisson measure 
on [0,oo) x [0, 1] x [— 7r,7r] with intensity measure dsda/3(9)dd. 

- Setting £ := E[V 2 ], for < e < 1, we consider the processes (V t ) t >o, (V t e ) t >o 
defined on (O, J", (J r t ) t > ,P) solutions to ([jD] and ([33]) with B, iV, W, W e defined 
previously, both starting from Vq. We also introduce the process (V^ e )t>o solution 
of the following S.D.E.: 



(5.9) 



(cos6- 1)V S £ _ - sm6W s -(a) N{dsdad0) 

Jo J\e\>e ' 

f f sin 6W S _ (a) N(dsdad9) 

Jo J\6\<e 



V!ds. 



By Proposition 3.1 and Proposition 3.3 C(V t ) and £(V/) are nothing but / ( and 



ft respectively. We set A e t =V t — V t e , At = V t e - V( and St (a) = W s (a) - W e s {a). 



Step 2: the aim is here to prove that 



(5.10) 



fi((A t e o ) 2 ) < d e e- Ce '° J ° e c ' s E a (S 2 s )ds + Ce 2 , 



where C depends only on £ and E[Vq ] and where c £ and d e are defined in (5.1). 
Making the difference between ( |5.9[ ) and (3.3), we get 

(5.11) 



Jo J\e\> 
t /.1 



(cos 9 - 1) Al_ - sin 9S e s _ (a) J N(dsdad9) - b e / A e s ds 
[ [ [ sm6W s -(a)N(dsdadO)- ^/2£b e B t 

JO JO J\8\<£ 



= H t 



A\_dX„ 



with (X t )t>o defined in (5.2) and with 

//, = - / / / sin edl_(a)N(dsdadO) 



JO J\8\> 
t r l 



f I I sindW s -(a)N(dsdad6) - y/2£b t B t 
Jo Jo J\e\<e 



i\e\< 

We do not write the dependence in e for H. According to Jacod [T^], = (L t 
D t )Z t , where Z t was defined in (5.4) and where 



(5- 



12) D t = - [ [ [ Z~}sm6W s -(a)N(dsdade)- y/2£b e [ Zj l dB s 

JO JO J\8\<e JO 
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and 
(5.13) 



JO J\8\> 



Z-}^-5 e s _(a)N(dsdad6). 
cos 9 



To verify this, it suffices to apply Ito's formula and observe that the process ((L t + 
D t )Z t )t>o satisfies the same S.D.E. than (A|) t > , i.e (L t + D t )Z t = H t + L{L S - + 
D s -)Z s -dX s . This S.D.E. has Lipschitz coefficients and thus has a unique solution. 
The processes (D t )t>o and (L t )t>o depend on e but we do not write this dependence. 
Hence 

(5.14) E[(A t e ) 2 ] = E[L t 2 Z 2 ] + E[D 2 Z 2 } + 2E[L t D t Z 2 }. 

- First, 

(5.15) E[L t D t Z 2 ] = Ex\E 2 (L t D t Z*)] - Ei[L t Z t 2 E 2 (A)] = 0, 

because for w = (wi,u; 2 ), we have (L t D t Z 2 )(u) — (L t Z^)(ux)D t (ujx, oj 2 ) and be- 
cause for uji fixed, E 2 [D t (wi, w 2 )] = 0. Indeed, recall that (L t )t>o and (Z t )t>o 
depend only on uj\ and that for wi fixed, N\Q\ <e (u]\, W2) is a Poisson measure while 
( J B t (wi,a; 2 )) t>0 is a Brownian motion on (f2 2 , J 72 , J" 2 , P 2 ), so that uj 2 )) t>0 
is a centered martingale (for uji fixed). 

- By Ito's formula, we have 



Z 2 L 2 t = -2b e I Z 2 L 2 ds 



o Jo J\e\>t 



[Z s _ + (cosfl - l)Z s _y ''(L s 
N(dsdade) 



smt 



cos ( 



zi Li 



2b e / Z 2 L 2 ds 



JO J\6\>t 



(cos 2 9 - 1)Z 2 S _L 2 S _ + sin 2 9{5l_{a)f 
- 2 cos 9 sin 6>Z s _i s _(5|_ (a)l N(dsdad9). 



Taking expectations and recalling (5.1 1, we get (use that cos 9 sin 0(3(9) d9 is odd) 

E(Z 2 L 2 )^-c e f E(Z 2 L 2 s )ds + d e f E a [(6l) 2 ]d s . 
Jo Jo 

Solving this differential equation, we find 

(5.16) E(Z 2 L 2 ) = d e e- c ' t f e c ^ s E a [{S e s ) 2 }ds. 

Jo 
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It remains to compute E(Z 2 Q D 2 g ). Recalling (5.12), we directly obtain 



E{Z 2 a D 2 o ) = E\(z tQ [ ° [ [ Z-}sm(-9)W s -(a)N{dsdad6) 



to 



Z t0 J Z- x ^fl£b t dB. 

V J0 JO J --it 
to 



(Z t0 Z~l ) ( Wl ) sin(-0) W s - (a)N$ <e (ua , dsdadO) 



(Z t0 Z- 1 )(u Jl )^/2£b:dB^(u J2 ) 



We thus obtain E(Z?D?) = Ei(W|(/x to , v to )), recall (5.6). We consider 



(5.17) 
with 



"to 



\8\<e 



sin 2 9/3{9)d9. 



Using the triangular inequality, we have 

wi k 1 > < ) < 2 ( wi (m^ 1 , c ) + wi (c , »c ) ) ■ 



By Corollary A. 2 and since E a (W 2 ) = £ for all s > by the energy conservation, 
we have (recall that u)\ is fixed) 



J t0 J^jZ^Z-yjuJr) sin 4 0E a (W?)/3(fl)d0rf s 
Jo° iiei<e(^*o^ 1 ) 2 ( w i) sin2 O£0(9)d9ds 



But E Q (W 4 ) < E(Vq 4 ) + 3£ 2 by Lemma A. 5 Furthermore, recalling (5.4), we have 
\Z ta Zj 1 \ = \e- b <( to - s) Y\ s<T .< t cos 9i\ < 1. Since finally sin 2 9 < 9 < e^on [-e,e], 
we easily deduce that for all lo\ fixed, 



C (E(U 4 ) + 3£ 2 ) 2 
£ 6 ' 



Finally, it obviously holds, recall (5.8) and (5.17), that for all oji fixed, 



WKCX 1 ) ^ jf" ( V 7 ^ - ^£ ^ sin 2 9~md9)\z- 1 Z t0 ) 2 (cj 1 )ds 



< £ 



2b f 



sm 



|9|<6 



(Z^Z^^ds. 



We used that (y/x - y/y) 2 < \x - y\. Recalling that \Zj x Z t <\ < e - b ^ to " s \ we 
easily get $1° (Zj 1 Z to f {ui)ds < ^. Furthermore, 26 e - L |<e sin 2 6>/3(0)d6» 

L, <e (2(1 - cos6>) - sin 2 9)fi{9)d9 < L |<e 9 4 (3(9)d9. Finally, it is easily checked 
that, since e € (0,1), & e = L, <e (l - cos9)P(9)d9 > \ f w<e 9 2 0(9)d9. Hence it 
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holds that for all oj\ fixed, 

f2( uJi cji 



< 2£e 2 . 



We conclude that Wf (m" S i^ 1 ) < Ce 2 (where C depends on £ and E[F 4 ]), whence 
(5.18) nzlD 2 tQ ] < Ce 2 . 



Gathering (5.14), (5.15), (5.16) and (5.18), we deduce (5.10). 



Step 3: in this step, we check that 



(5.19) 



E[(A<) 2 ] < -e 2 



for alH > 0. We first observe that ( |3.1[ ) can be rewritten as 

ft rl 



V t =V + 



JO J\9\> 
t r l 



JO J\9\<e 



(cos6- l)V a - - sm6W s -(a) N(dsdad6) 



(cos6- 1)V S - - sm6W 8 -(a) N(dsdad0) - b t / V s ds 



Hence, making the difference with (5.9), we find 

ft pi 

(5.20) 



A £ t = f [ j (cosd -l)A e s _N(dsdad9) 
Jo Jo J\e\>t 

+ [ [ [ (cos9- l)V s -N(dsdad6) -b e [ A f s ds. 

JO Jo J\9\<e Jo 



'0 JO J\9\<e 

Applying Ito's formula, we get 



o Jo J\9\< 
t r i 



(A e t ) 2 = -2b e f (Al) 2 ds+ f f f (cos 2 6»- \)(A £ s _) 2 N(dsdad9) 

Jo Jo Jo J\6\>£ 

(a € s _ + (cos9- 1)V;_) 2 - (A^_) 2 N(dsdade) 

Al + (cose-l)V s ) 2 -(Al) 2 
- 2A^(cos6» - l)v}p(6)dsdad0 
(2b e [ (Al) 2 ds+ f f f sm 2 e(A f s _) 2 N(dsdade)) +M t 

^ Jo Jo Jo J\9\>e ' 

V 2 (l- cos 0) 2 l3(e)d6ds, 



JO J\9\<e 



10 J\9\<<l 

where (M t ) t > is a centered matingale. Taking expectations, this yields, recalling 
(5.1) and that E[V S 2 ] = £ for all s > by the energy conservation, 

rt 



Thus 



E[(A?) 2 ] = - Ce / E[(Al) 2 ]ds + le £t. 
Jo 

-r.t\ , 



™r,A f sQi le£ -r ts - 
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But 7e < (because for |0| < e, (1 - cos 9 ) < 9 2 /2 < e 2 /2) and c £ > 2b e . We 
deduce that 7 c /c e < e 2 /4 and finally get (5.19). 

Step 4: we now check that E[AjA^] = for all t > 0. Applying Ito's formula, 
using (5.201 and (5.11), we have 



AidB, - 2b e 



AlAlds 



J ^|9|>e 



(cos0A|_)(cos0A|_ - sm0<^_(a)) 



- A!_A! 



N(dsdad9) 



III \(A e s _ + {co S e-l)V s _)(Al_- S meW s -(a) 

JO JO J\6\<e LV 7 



A!_A! 



N{dsdadB) 



III \ ( Al + (cos 6 -l)V s )(Al- sin 9W s (a) 

Jo Jo J\6\<e LV ' 



- A\A\ - A^(cos0 - 1)V S + A £ s sm9W s (a) f3(9)dsdad9. 
Taking expectation and using that j3 is even, we get 



EfAjAf] = -2b c / E[A«A*]ds 



(cos 2 9 - l)E[A e s A e s ]f3(9)d9ds. 



o J\e\>t 



So the function t ^ E[A|A|] solves the O.D.E. y' = -c e y, see (|5_lJ). Since y(0) = 0, 
we easily conclude. 



Step 5: conclusion. Using Steps 2, 3 and 4, we find that 

mVt - V: a ) 2 ] =E[(A^) 2 ] +E[(At o ) 2 } + 2E[A^ o Ay 



< d e e- c ' to 



s E a [(5l) 2 ]ds + Ke 2 , 



where K depends only on £ and E[V 4 }. We set u(t) = W|(/t>/t e ) = E Q [(5|) 2 ] by 
@. Since C(V t0 ) = f t „ and C(V£) = f? , we have u(t ) < E[(V to - VQ 2 ]. Since 
to > is arbitrary, we get, for all f > 0, 

u(t) < d e e~ c ^ I e c ' s u{s)ds + Ke 2 =: v(t). 



Consequently, 

v'(t) = -c e (v(t) - Ke 2 ^j + d e u(t) 

< -c 6 (v(t) - Ke 2 ^j + d e v{t) 

< (d € - c e )v(t)+c e Ke 2 . 

We first observe that d e < c e , so that v(t) < v(0) + Kc e e 2 t < Ke 2 (l + c e )(l + t) < 
Ce 2 (l + 1), because c e < /* 2 /3(0)d0, see (pul 
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We can also obtain a uniform in time bound. Recall that v'(t) < (d e — c e )v(t) - 
c e Ke 2 = — 2b t v(t)+c e Ke 2 . We observe in fact that v'(t) < as soon as v(t) > C| /v ' 



26 e 



Since v(Q) = Ke 2 < &^pS we classically deduce that v(t) < < C e v for all 

t > 0. 



So we have W 2 2 (/t,/t) = u(t) < v(t) < Cmm (e 2 (l + t), f-\ for all t > 0. To 
complete the proof of Theorem 2.5 it suffices to observe that 46 e > Jj 9 | <e 9 2 f3(9)d9 
for any ee (0,1). □ 



6. Convergence of the particle system 

In this section, we prove the results about the approximation of the solution of 
the Kac equation by a system of particles. Let thus fo G V^M.) and let /3 be a cross 
section satisfying (2.1). We fix e £ (0, 1), and we consider the solutions (/t)t>o and 
(/t)t>0 to (1.1 1 and (1.51 respectively, both starting from f . 



In the first part, we will rewrite the system of particles (2.11 1 in a suitable way 



and in the second part, we will introduce a system of i.i.d. particles with law 
(ft )t>o- Using these systems of particles, we will be able to prove Theorem 2.7 and 
its corollary. We will end this section with the proof of Proposition 2.9 and with 
an extension about the Wasserstein distance W 1 for 7 E (1, 2). 

We recall a usefull result (see e.g. Villani [25, Remark 2.19 (in)])- 

Proposition 6.1. If /i and v art two probability measures on K, for 7 > 1, we have 
W^Jlfiji 1 ) = Jq 1 (F f f 1 (a) — i< 1 I 7 1 (a)) 7 G?a where F fJi (x) = fi{(— 00,2;]) and F v (x) = 
^((—00, xfj . 

6.1. Another way to write system ( |2.11[ ). We fix an integer n and we consider: 

• a family of i.i.d. random variables (Vq)^^ n \ with law f , 

• a family of i.i.d. Poisson measures (iV 4 ) i6 { li on [0,oo) x [0,1] x [— tt,it] 
with intensity measure dtdaf3(9)d9, 

• a family of i.i.d. Brownian motions (-Bt)t>o,ie{i,...,n} ■ 

For e € (0, 1), we consider (V t l '™' e ) t > , ie{i,...,«} solution of the following system of 
S.D.E.s: for i = l,...,n, for all t > 0, 



vr 



-V 1 



JQ J\8\>e 
t 



(cos0- 1) VI 



sin6» F, 



N l (d6dads) 



b £ f V s hn ' e ds + y/2£b~eBl 
Jo 



where F^ = i E?lv?--<x' £ = Jk^/oW and b t = Jj fl|<£ (l - co S 6)[3(8)dd 



This particle system is identical (in law) to the one introduced in (2.11). In- 



deed, it suffices to note that given (^^"'"Oie-fi,...,™}) the law of (F™^) 1 (a) (with 
a uniformly distributed on [0,1]) is the same as that of V/l"' 6 (with j uniformly 
distributed in {1, ...,n}). 
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6.2. A system of i.i.d. particles. For i 6 {1, ...,n} and e > 0, we consider the 
process (V"/' e )t>o solution of the following S.D.E. (with the same random objects 
Vq, N l and B % as previously), 



JO J\8\>e 
t 



(cos9- 1)V^ - sinflfeO 1 (a) N\dsdad6) 



Jo 



where F?(x) = £ 00 ft (dv). 

For each s > 0, it holds that CauFg) -1 )) = / e . Hence we can apply Proposition 
with W s £ = (F^)- 1 and deduce that for each i € {1, n}, each t > 0, £(V^' e ) = 



3.3 



/ t £ . Furthermore, the processes (V^)t>o are obviously i.i.d. (for i = 1, ...,ra). 



6.3. Proof of Theorem 2.7 , We start with the following result. 

Proposition 6.2. We set /z™' e = i Ei $v i,n,€ an< ^ Pt = n Si <V ie ■ Then for 
any T > 0, 



sup E W^(/ t e ,Mr) < C(l + T) sup E ^(/ t e , 

[0,T] L J [0,T] L 

where C depends only on Jq 9 2 (3(9)d6 . 

Proof. To lighten notation, we set V/ = V'/'™' e and VI = V^' e for the whole proof. 
By the triangular inequality, we have W^(/t , < ^(/t 6 , M™' 6 ) + ^(/i™' 6 , M™' £ )- 
Hence, by squaring and taking expectations 



(6.1) 



E 



-2E 



w 2 (/ t %Mr)w 2 (Mr,Mr e ) 



Using the fact that W| ( ~ Ei $e< > £ Ei O ^ S El \ x i ~ Vi ? , we have 



E 



wnpr, & £ j\ < e(- e - ^i 2 ) = E (i^ - y t 



i=X 



rl|2 



We set A t = V^ 1 - V^ 1 . It holds that 



(cos0 - 1)A S _ - sin0<5 s _(a) ^(dsdadO) -b e A s ds 



'o Jo J|e|>c 

where J* (a) = (F" ,e ) (a) — (Ff) (a). Applying Ito's formula, we get 



J0 J|9|>e 



(A s + (cos (9- 1)A S - sin65 s (a)) 2 - A 2 S N 1 (dsdadO) 



2b f / Aids. 
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Taking expectations and using Proposition |6.1[ we get, with c £ and d e defined in 
(PI 



«(*) :=E(A?) 



f E{A 2 s )ds + d e E^jf <5 2 (a)da)ds 



( S )ds + d £ / E[W£(/. e ,/# e )]d 



Jo 



the last equality being obtained by solving the differential equation satisfied by u. 



If we set u(t) := E Wf (j?, /i"'^ 



, s n := sup [0 T] E 
0,T], 



and if we return 



to (6.1 1, we thus find, for all t € 

u (t) < s n + v (t ) + 2E [\V 2 {ft , ) W 2 ( ^ , M r ) 

by the Cauchy-Schwarz inequality. We thus have 

«(t) = d £ e- c «* / e c » s M(s)ds 
Jo 
ft 

<d f e~ c ' t 1 - c ' s 



< s n + d e e 



(s n + v(s) + 2y/s^y/v(sfj ds 
e C€S (v(s) + 2^/s^^/v{s)^ds =: w(t). 



We used that d £ e Ce * f*e c ' s ds = ^(1 - e Cet ) < 1, since c? £ < c £ , recall (5.1). 
Differentiating w, we get 

u/(t) = -Ce(w(t) - 8 n ) + d e (v(t) + 2y/^y/v(t)) 

< c e s n + w(t)(—c e + d e ) + 2d ey /s^^ w(t) 

< c e s n + 2d ey r s^^w(t) 

< as n + 2a y /fh n y / w(t), 

where a = 9 2 (3{e)d6. 

Putting x(t) = w(t)/s n , we deduce that x(0) = 1 and x'(t) < a(l + 2yJ x(t)) < 
2a(l + y/x{t)) < Aa-y/l + x(t), whence y/l + x(t) - y/l + x(0) < 2at, which gives 
x(t) < C(l + t) 2 , and so w(t) < C(l + t) 2 s n . 

To summarize, we have v(t) < w(t) < Cs n (l + 1) 2 and, for all t <E [0, T], 



E 



W 2 2 {fl,$' e ) = «(*) < + «(*) + 2vW«(*) 



< s„ + Cs„(l + t) 2 + 2y/^VCv^{l + t) 

< Cs n (l + t) 2 

= csu P E[v^ 2 a^r)l(i+;) 2 - 

[0,T] J 



This concludes the proof. 



□ 
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Theorem |2 . 7| follows almost immediately. 



Proof of Theorem \2.7\ For each t > 0, we consider an i.i.d. sequence 
{Vt)i £ {i,...,n} with law f t such that for each i, E[(V?-^f' e ) 2 ] = W£{f t ,f[). Then we 
set /j," = Using the triangular inequality, Theorem 2.5 and Proposition 

6.21 we have 



sup E 

[0,T] 



(6.2) 



Wi{f u ^) < 2 sup Witft, ft) + 2 sup E WiUt, ^) 

J [0,T] [0,T] L 

< C(l + T)e 2 + C(l + Tf sup E[w2(/ t e ,#' 6 ) 

[0,T] L 



We use again the triangular inequality to obtain 



E 



But E 



5J we get 
(6.3) 



(tf,tf' e )} < £E?=l E (l^-i?'T) =Wi(f u ft). So using Theorem 



sup E 

[0,T] 



^l(/t e , /V )1 < C((l + T)e 2 + S upEfw|(/ t ,M?) 

J V [0,T] L 



Inserting (6.3 1 in (6.2), we obtain 



supE Wi(f t ,tf> e ) <C(1 + Tf e 2 + su P E W|(/ t ,#) 

J V [0,T] L 



[0,T] 

which concludes the proof. 



□ 



Finally, we give the proof of Corollary |2.8 



Proof of Corollary |2.8| It suffices to apply Theorem |2.7[ Lemma |A.4| with 
7 = 2 and q = p — 7, and Lemma |A.5| □ 



6.4. Other Wasserstein distances. The first part of the following result is Propo- 
and in the second part, we give some estimates about E[W^(/t, /u"' e )] for 



sition 



2.9 



Proposition 6.3. Adopt the same notation as in Theorem pT 
(i) If we assume that fo € V4,(M.) and if 8(5(9)d8 < 00, then 



B apE[W 1 (f t ,£' e )]<C T (e+- 

[0,T] V V 

where Ct depends only on T, f and /?. 

(H) V fo € PpQR) /or some even p > 4 and if ft 6"* P{6)d6 < 00 for 
7 € (1,2), tfterc 



mipEp^C/t,^)] < C T (e 7 + 

[0,T] v n 2 P- 2 



where Ct depends only on T, f Q , j3, p and 7. 
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Proof. Let 7 £ [1, 2) be fixed. We assume that J* 9^/3(9)d9 < 00 and we set 
V t z — V/'™' e and V t l = U/' e for the whole proof to lighten notation. 
Step 1: we first prove that 

I 7 , e (jc):= / [|cos6'-a;sin6ip-l]/3(6i)^<C(l + |x| 7 ). 
J\e\>e 

To this end, we set J 7 , e (a;) = /| | >e [|1 — £#| 7 — l]/3(0)d#. Using the inequality 
| 7 _ 57 1 < q _ foKaT- 1 + feT- 1 ), we get 

|/ 7£ (a;) - J 7 e(x)\ <C I \(cos9-l)-x{sm6-6)\ 
J\e\>e 

(I cos (9 - xsindp- 1 + |1 - atflT- 1 )^)^ 



'/ 2 (l + |a;|)(l + |sp r - 1 )0(0)d0 



< c 

/|0|>e 

^^(l + lx^l + lxp- 1 ) 
<C7(l + |x| 7 ). 

Using the fact that /3 is even, we can write 



with 



and 



J-yA*)= I [\l - x8p - 1 + 7 x9]f3(9)d6 = J*(x) + J 2 (x) 

J\6\>e 

J^ e {x)= f [\l-x9\^ -l + -fx9]p{9)d9 

J\9\>e,\x6\<l/2 



J 2 e (x) = ( [\l-x9p -l]p(8)d9. 

J\e\>e,\x6\>l/2 

By Taylor's formula, we get, observing that \x 2 9 2 \ < |x| 7 |6»| 7 if \x9\ < 1/2, 
J 7e (x)<C f x 2 9 2 (3(8)d9 <C\x\~< f \9\i/3(9)d9. 

J\6\>e,\x0\<l/2 J-tt 

Next, since 1x6*1 > 1/2 implies 1 + \x9\~> < (1 + 2 7 )|a;(9| 7 , 

e (x) <C [ [1 + \x9\~<](5{9)d9 < C|x| 7 f \9\~<f3{9)d9. 

J\8\>e,\x6\>l/2 J-tt 

We thus have J 7 , e (a;) < C\x\ J and hence I lt e{ x ) — ^(1 + \%\ J )- 



Step 2: using Step 1, we now prove 

W]{ff e ,tf e j\ <Ce ct f E [W 7 (/|, ds, 

r W7(^ ,£ ,/i? ,e )l <E(|A t | 7 ), where A t = V t 1 -V t 1 . It holds 



E 

for all t > 0. We have E 



A t =/ / j [(cos 6* — 1)A S _ — sin 95 s - (a)] N 1 (dsdad9) — b e f A s ds, 
Jo Jo J\e\>e L - 1 Jo 



10 Jo J\0\>e 
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where 8t{a) = (i* 1 ™' 6 ) (a) — (-F t e ) (a). By Ito's formula, 

[| cos6»A s - sirL0S s (a)p - \& s \"t]p(8)dsdad6 

7 6 £ / E(|A a p)ds 
Jo 











E 




= E 


[/77 








-Jo Jo J\ 



< E 



|0|>e 



[|cos0-sin/ s(a)| 



A. 



i\p(6)dsdad6 



Using Step 1 and then Proposition [Ol we get 



E 



|A,P 



< CE 



I A, P I 



] dads 



< C / E[|A fl | 7 ]ds + C 
We conclude by Gronwall's lemma. 



E 



ds . 



Step 3: using very similar arguments as in the proof of Theorem 2.7 
m£tort~W y (f t , ft) < W 2 (f t , / t e ), we get 



and observ- 



E 



E 



for all t > 0, where //" is the empirical measure of a sequence of i.i.d. random 
variables with law ft- 

Step 4: the aim of this step is to prove that 



(6.4) 



sup E WJ{f t , tf> e ) <C T eH sup E W](f t , tf) 



[0,T] 



[0,T] 



Using the triangular inequality, Theorem |2.5[ Step 2 and Step 3, we have 



< 



€< +E 



Ct 



E 



ds 



Using again the triangular inequality and Theorem |2.5[ we get 



E 



We thus have 



w^ft^T, 



< C E 



E 



w^(f s ^r) 



< C(l+T)^ 2 eHCE[w(/ s ,^ e ) 



E 



w](f t ^r) 



< C T (e 



E 



W](f t ,tf) 



E 



w]{fs,nV) 



ds), 



and we conclude with the help of Gronwall's lemma. 



Step 5: we can now prove (i). Since /o € ^(K), Lemma A. 5 implies that 
su Pt>o f R vi ft(dv) < oo, whence sup t > E Wi(f t ,(4) < ^= by Lemma 
sorting this in (16. 41), we easily conclude. 



A.3 



In- 
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sup t > / R v p f t (dv) < oo, whence sup t > E (ft, fit) < P - T by Lemma 

J n 2p-2 



Step 6: we finally prove (ii). Since fo € 'Pp(R), Lemma A. 5 implies that 

7 



Inserting this in (6.4), we easily conclude 



A.4 



7. Numerical results 

We consider here the cross section (3(9) = with < v < 2. Let /o be a 

probability measure admitting a moment of order 4. We fix an integer n, a small 
parameter e > 0, and we take the same notation as in Section |2.4| 

We simulate two systems of particles : the system (V/'™' £ ) t >o.ie{i,...,n} described 
in Section |2.4| (system with diffusion) and the following system without diffusion: 
for i £ {1, n} and t > 0, 

K n,t =V r i + /* / / [(cos^ - l)f;i"' e - sin^/l"' e l^'™(d S dMj). 
Jo Jo J\e\>e L 

The algorithm is the following (we write in italic the parts which only concern 
the system with diffusion). 

We set t = 0, and for i = 1, n, we simulate V(i) ~ /o and set T up (i) = 0. 
While t < T /i,lai (where T /i,lai is the time that we want to reach), we 
simulate an exponential random variable T with parameter n Jigi >e P(9)d6 
and we put t = t + T. We choose randomly two integers i and j in {1, n}. 
For o«r system with diffusion, we update the particles i and j by setting 

V{i) = V(i) cxp(-6 e (i - T^(i))) + G(i), 

and 

y(i) = y(j) exp(-6 £ (i - T^(j))) + G(j), 

where G(i) (resp. G(j)) has a centered Gaussian law with variance 1 — 
exp(-2b e (t-T w P(i))) (resp. l-exp(~2b e (t-T u P(j)))), whereb, is defined 



m (5J\ ), and we set T u P(i) = T up (j) = t. 

Next, we simulate a random variable with density /3 e /||/3 e ||i, where 
(3 e (9) = /3(0)li0i> e , and then, for the two systems, we put V(i) = cos<dV(i) — 
sin6F(i). 

Only for the system with diffusion, we update all particles with a Gaussian 
term: for i 6 {l,...,rt}, 

V(i) = V(i) exp(-b e (t - T up (i))) + G(i), 

where G(i) has a centered Gaussian law with variance 1 — exp(— 2b e (t — 
T up (i))). 

For our simulation, we take T final = 0.1. Our initial data is f = + 5 x )/2. 
The goal here is to see what system is more efficient. For this, we need a reference 
curve. We obtain it by simulating n = 10 7 particles with e = 0.03, and by using a 
smoothing procedure. 

We see that the system with diffusion term is much more efficient when v is close 
to 2. For v smaller, the difference not clear. 
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Figure 1. v = 0.5. Left graphic: system with diffusion, n = f0 4 , 
e = O.f. Right graphic: system without diffusion, n = 2.f0 4 , 
e = O.f . Both simulations need approximately 0.05s. 




Figure 2. v = 0.5. Left graphic: system with diffusion, n = 10 4 , 
e = 0.1. Right graphic: system without diffusion, n = f0 4 , e = 
0.02. Both simulations need approximately 0.05s. 




Figure 3. v = 1.5. Left graphic: system with diffusion, n = 10 4 , 
e = O.f. Right graphic: system without diffusion, n = 2.10 4 , 
e = O.f. Both simulations need approximately 0.f4s. 

Appendix A. Appendix 

A.f . Wasserstein distance between a Poisson integral and a Gaussian law. 

We start with a result of Rio [TSJ Theorem 4.1], which gives some very precise rate 
of convergence for the standard central limit theorem in Wasserstein distance. 
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Figure 4. v = 1.5. Left graphic: system with diffusion, n = 10 4 , 
e = 0.1. Right graphic: system without diffusion, n = 10 4 , e 
O.Ofi. Both simulations need approximately 0.14s. 



Theorem A.l. There exists a constant Co such that for any positive integer n, 
for any sequence (Xi)i>o of real independant centered random variables in L , 



where r\ n = C(v n S n ), S n = Yh=i X *> Vn = Var(S n ). 

Using this result, we can estimate the Wasserstein distance between a compen- 
sated Poisson integral and a centered Gaussian law with the same variance. The 
following result is very close to Corollary 6]. 

Corollary A. 2. If E is a Polish space endowed with a non-negative a -finite mea- 
sure v, if N is a Poisson measure on [0, T] x E with intensity measure dtv(dz) 
and if H : [0, T] x E ^ R is a deterministic function such that J Q J E (H 2 (s 7 z) + 
i? 4 (s, z))v(dz)ds < oo, then setting 

X t = [ [ H(s,z)N(ds,dz), q t = [ [ H 2 (s, z)v{dz)ds, 

JQ JE JO JE 

we have 

o, , L H i (s. z)v(dz)ds 

It 



where Cq is a universal constant (the same as in Theorem A.l). 
Proof. For n > 1, i G {1, ■•■,«}, we set 

fit/% 



X? = Jn 



/ / H(s,z)N(ds,dz) and S„ = Vl,' 

J(i-l)t/n JE 



We have X t = The random variables X" are independent, centered, 

pit / n t> n 

E[{X?) 2 ]=n H 2 {s,z)v{dz)ds, v n = Var(S n ) = V E[(Xf) 2 ] = 

J(i-l)t/nJE i=1 

It classically holds that 

rit/n /> pit/n /> 2 



nq t . 



E[(Xf) 4 ] = n 2 / / H 4 (s,z)v(dz)ds + 3n 2 ( / H 2 {s,z)v{dz)ds 

J(i-l)t/nJE \J(i-l)t/nJE 
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Hence 



« ft r 

^E[(A™) 4 ] =n 2 / H i (s,z)iy(dz)ds 

™ /-it/n 2 

+ 3n 2 V( / / H 2 (s,z)v(dz)ds) 

j-[ K J{i-l)t/nJE ' 



By Theorem [Al| 



(C(X t ),M(0,q t )) = W 2 (c(-^=S n ),M(0,q t ) 
= q t W 2 (c(-^S n ),Af (0,1) 



< 



J*J E H 4 (s,z)v(dz)ds 

.On / / 



C / r lt / n r s 2 

9t i=l ^J(i-l)t/n JE ' 



Setting F(f) = J Q |/(s)|ds with f(s) = J E H 2 (s, z)v(dz) and observing that F is 
continuous (and so uniformly continuous on [0, T] for all T > 0), we obtain that 

S"=i f /(i— "){/n /b ^ 2 ( s ' z)v(dz)ds^j — > when n — > +oo. Since the last formula 
holds for all n > 1, we easily conclude. □ 

A. 2. Rate of convergence of empirical measures. We first give a classical 
result about the Wasserstein distance W±. 

Lemma A. 3. Let [i be a probability measure in V4QH) . We consider n i.i.d. random 
variables (Xi) ie {i ....,„} with law /1 and we set fi n = rX^iLi^*' Then there exists 
a constant C depending only on J R x \i(dx) such that 

L J y/n 

Proof. If we set F(x) = n((— 00, a;]) and F n (x) — ^Yli lx f <n we bave (see 
Villani Hi p 75]) 



e(w x (m,M„)) =E( / |F(a;) - F n (x)\d3 

/oo ^ n 

\F(x)-~^2l Xi < x \di 
-00 ri • 1 



If Y ~ B(n,p), e(|£ -p|) < Je((^-p)2) = ^Vp(l-p). Hence, since for 
each x, £)" =1 lxj<x ~ #(n, F(af)), 

/ \ 1 f°° 

^[W^,^)) < — y/F(x)(l-F(x))dx. 
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But A := j„x fi(dx) < oo implies that for x > 1, (1 — F(x)) = /Lt([x, +oo)) < A/x* 
and for x < -1, = n((-oo,x]) < A/x A , so that y/F{x)(l - F(x))dx < 

oo. □ 



We now deduce similar estimates for other Wasserstein distances. 



Lemma A. 4. Let n be a probability measure admitting a moment of order (7 + 7, 
with 7 > 1 and q > 0. We consider n i.i.d. random variables (-X"»)jg{i n } with 
law 11 and we set fi n = ^ • There exists a constant C depending on 7, q 

and on the moment of fi of order (7 + 7 such that 



E 



W7(/i,Mn) 



< 



c 



^2(9+7-1) 



Proof. Let us denote by (f2, J 7 , P) the probability space on which X\, X n are 
defined. For a fixed w £ fi, we consider two random variables X and Y u defined 
on the probability space ([0, 1],2?([0, l]),da) with C a (X) = fx and C a {Y u ) — n n {w) 
such that Wt(p, n n (u))) = E a (\X - Y^l). Then we have, for any A > 0, 

Wjfa m»M) < E Q (|X - r-H < A^W^n, » n (oj))+E a (\X - Y"yt lx _ Y „ l>A ). 

We observe that 



E a (\x -Y u \ q+ "f) 



E a (\x-Y»\n {x _y„ >A ) < 

But, setting m p (/i) = J R 

E a (|X - FT +7 ) < C(E Q (|X|^) + E Q (|y<T +7 )) = C(m ?+7 (/i) + m g+7 (,i„H). 



One easily checks that E m g + T (/x n ) 



(/i) . Using Lemma 



A.3 



we finally get 



E 



< C 



^ 7 1 ™g+ 7 0) 



Choosing A = n^+T- 1 ) completes the proof. 



□ 



A.3. Moments of a solution to ( fTTT| ). In many places o f the proof, we need to 
upperbound J R v 4 f t (dv) for any t > where {ft)t>o solves (1.1 1. We also need to 
upperbound higher moments. 



Lemma A. 5. For f € T-^R), consider the unique solution (ft)t>o t° {1-1). For 
any t > 0, we have 



v 4 f t {dv) < / « 4 /oW+3 / v 2 f (dv) 



If fo G V P (M.) with p even, then there exists a constant C depending on p, fj ana 
on J R v p fo{dv) such that for any t > 0, 



v p f t (dv) < C. 
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Proof. We only treat the case p = 4, see Truesdell |22| an d Desvillettes [5] for 
the general case. If we take (p(v) — v 4 , we find, recalling (2.2 1 and using that j3 is 
even, 



K%(v,v*)= / (cos 4 9 - l)v 4 + sin 4 9vt + 6cos 2 9sm 2 9v 2 vi P(9)d9. 

J —TT 

Setting m&(/i) = f R v k n(dv) for /i a probability measure on K and k G N, we thus 
get, using 

"T- 4 (/t) = m 4 (/ ) 



/ f [- (1- cos 4 0- sin 4 9)m 4 (f s ) 

JO J-TT 



+ 6 cos 2 9 sin 2 9m 2 2 (f s )]f3(9)d9ds. 

Recalling that ni2(f s ) = ru^ifo) ='■ £ for any s > 0, observing that (cos 2 - 
sin 2 9) 2 = 1, which gives 2 cos 2 sin 2 9 = 1 — cos 4 9 — sin 4 9 and setting c = /^(l • 
cos 4 - sin 4 9)f3(9)d9, we have 

TO 4 (/t) = m 4 (/o) - c / m 4 (f s )ds + 3c£ 2 t, 
Jo 

whence 



m 4 (/ t ) = (m 4 (/ ) - 3£ 2 ) exp(-ci) + 3£ 2 < m 4 (/ ) + 3£ 2 



as desired. 



□ 



sup I (h 2 + 2vh)q(t,v,dh) <C T (l + v 2 ), 

[0,T] 



A. 4. Well-posedness for a P.D.E. To conclude this paper, we state the following 
result. 

Proposition A. 6. For t > and (v, v*) G K 2 , we consider two finite non-negative 
measures q(t,v, dh) and r(t,v,v*,dh) onM. such that A q := sup t v q(t, v, K) < oo, 
A r := sup t v Vt r(t, v, v*, M) < oo and /or aZZ T > 0, all (v, v*) G K 2 , 

(A.l) 

and 

(A.2) sup / (h 2 + 2vh)r(t,v,v*,dh) < C T (1 

[o.T] Jr 

Let also a > and b G M &e fixed. Then, for any /o G 7^2 (K) > ^ftere exists a unique 
(ft)t>o G L^ c ([0,oo),P 2 (]R)) smc/i ttai/or aZZ <p G C b 2 (M), a// 1 > 0, 

^ f v(v)f t (dv) =a [ <p"(v)f t (dv) + b [ v<p'(v)f t (dv) 
(A.3) +/ / \<p(v + h)- <p(vj\ q(t,v,dh)f t (dv) 

ip(v + h) - <^(u) r(t,v,v*,dh)f t {dv)ft(dv t ). 



+ v 2 + v 2 ), 



Proof. We denote by A^(M) the set of finite signed measures on K. If \x G A^(M), 
we set \h\tv — su P ¥ >gL°°,|| ¥ )|| 00 <i Jr <p( v ) ^{dy) ■ Using the Lusin Theorem (see e.g. 
(U Theorem 9.11]), we have \h\tv = su PyeCi,,|M|oo<i Jr L P( v )l I (^ v )- We also have 
ImItv = Jr ImK^) where = /i+ + /i_ and if /i has a density / with respect to 
the Lebesgue measure, \/j,\tv = J R \f(v)\dv. 
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Preliminaries. For e > 0, we set G e (v) 



/2we 



-e 2« 



Let fie M(R). We claim 



that lim e ^o Im * G £ \tv — |m|tv- Observe that this is not obvious, since it does not 
hold true, generally, that lim e _ >0 \v*G € — fi\rv — (choose e.g. fi = 5 ). First, we 
have 



\fi*G t \Tv = / I A* * G e (v)\dv = / / G e (v — w)n(d 

JR il JR 



dv 



< 



G e (v — w)dv\fi(dw)\ = \fj,\ 



TV- 



Next, let a > 0. There exists a function ip e Cb with ||<p||oo < 1 such that 
J R (p(v)fi(dv) > \fi\rv — a - We have, since \x * G e clearly converges weakly (in the 
sence of measures) to fi, 



\h*G c \tv> I <p{v)(/j,*G e )(v)dv 



(p(v)fj,(dv) > \h\tv - a. 



Making a tend to zero, we get limmf e ^ Im * G e \rv > ImItv- 

Uniqueness. We consider two solutions (ft)t>o an d (<7t)t>o ; with f = g and for 
t > we set fj, t = f t — g t . For any <^ e C^(M), any t > 



d 
eft 



(p(v)fi t (dv) =a / ip"(v)n t (dv) + b / vip' (v) fi t (dv) 

JR JR 

+ <p(v + h) — <p(v) q(t,v,dh)nt(dv) 

iiii L - 1 

+ ft) - r(t,v,v^dh)[f t (dv)ft(dv^) 

- g t (dv)g t (dv*)]. 



+ 



We first observe that f t (dv)f t (dv*)-g t (dv)g t (dv*) = ft(dv)fj, t (dv*)+g t (dv*)(j, t (dv). 
We have 



d t {p. t * G e )(v) 



dt 



G e (v - w)fit(dw) 
a / G"(v — w)fi t {dw) — b / wG' e (v — w) /i t (dw) 

JR. JR 

/ / [G e (v — w — h) — G e (v — w)]fit(dw)q(t, w, dh) 

JR JR 

f f f [G e (v - w - ft) - G e (v - w)}[f t (dw)fi t (dw*) 

JR. JR JR. 

+ g t (dw*)fj, t {dw)]r(t, w, w*,dh). 



For t) > 0, we consider a function T v of class C 2 such that for any (\x\ — i]) + < 

T v (x) < \x\, r;j(x) > and Hr^l^ < l. We also assume that T Vl > T m if ^ < rj 2 - 
Observing that J R G"(v — w)/j, t (dw) = (fi t * G e )"(v), we have 



d_ 
dt 



f r n ({iH*G e )(v))dv= f r , v {{f H *G e )(v))d t (iH*G e ){v)dv 

JR. JR 

= A t + B t + C t + D t , 
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where 

A t = a [ T'^int * G e )(v))(fH * G e )"(v)dv, 

JR. 

B t = -b f f r 7 ((/^ * G e )(v))wG' e (v - w)iH(dw)dv, 

JR JR 



Ct= / / r^((/i t * G e )(v))[G e (v — w — h) — G e (v — w)]/j, t (dw)q(t, w, dh)dv, 

JR JR JR 

and 

D t = [ [ [ [ T' v ((p t *G e )(v)){G e (v~w-h)~GAv-w)}{f t (dw)fi t (dw*) 

JR JR JR JR 

+ gt(dw«)nt(dw)]r(t, w, u>*, dh)dv. 
Using an integration by parts and recalling that T^'(x) > for any x e M, we have 

A t = -a [ r;'(( Mt * G e )(v))((jH * G € )'(v)) 2 dv < 0. 

JR 

First writing w = v + w — v and then using an integration by parts (observe that 
J R G' e (v - w)m{dw) = (fi t * G e )'(v)), we have 



B t <-b / r;(( Mt * G e ){v))vG' e {v - w)nt(dw)dv 

JR JR 

+ l|r;i|oc|6| / / \w - v\\G' e (v - w)\\^ t (dw)\dv 

JR JR 

<b f r v {{to*G e )(v))dv + \b\ ( \v\\G' e (v)\dv f \tH\(dw) 

JR JR JR 

< C( J \(txt * G e )(v)\dv + I^Itv) < C\ih\tv- 

We used the preliminaries and the fact that J R \v\ \G' e (v)\dv < C. Using next that 
J R G e (v ~ w — h)dv = J R G e (v — w)dv = 1, we have 

C t + A<2||r;||ooA, / |A* t |(d«;) + 2||r / || 00 A r ( f \nt\(dw.) + f \vt\(dw)) 

JR K JR JR ' 

< C\fit\ T V- 

We thus get 

d 



dt 



T v (nt*G e {v))dv<C\Lit\ 



TV- 



Using the monotone convergence Theorem (recall that T v (x) increases to |x| as rj 
decreases to 0) and recalling that fi = 0, we have 

\lH * G t \ T v = lim / T v ((fi t * G e )(v))dv 

< lim / ^((no *G e )(v))dv + C / \fi s \ T vds 

V-tOjR JO 

<G \/j, s \ T vds. 
Jo 
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Making e tend to and using the preliminaries, we get, 



I I tv < C \ \fj, s \ TV ds, 
Jo 



and we deduce that |/*t|TV — by Gronwall's lemma. 

Existence. For (Qt)t>o € Lf£ c ([0, oo), ^(K)) given, we consider the following 
linear P.D.E. with unknown {g?) t >o- for all tp e C^(M), all t > 0, 



(A.4) 



<p(v)g?(dv) =a I V "{v)g?{dv) + b / vy'(v)g?{dv) 

JR 

(p(v + h) - <p(v) q{t,v,dh)gf(dv) 

ip(v + h) — ip(v) r(t, v, v*,dh)gf(dv)Qt{dv*). 



For t > and € R 2 , we consider the following probability measures 



Vt,v(dh) 



q(t, v, dh) 



A n 



and 



r(t, v, v*,dh) 
A r 



A r 



So(dh), 



and we set F? iV (x) := ^((-oo, a;]) and ^ (x) := ^^((-00,2;]). We also 
set H q (t,v,u) := (F J ? 1 ,) _1 (u), H r (t,v,v«u) := (F^ and we consider the 

following S.D.E. 



r;=I -/ / H q {s 7 V s ^,u)N q (dsdu)+ [ [ [ H r (s,V s ^,v*,u)N r (dsdv*du) 



^0 



JUO 



(A.5) 



6 / V s ds + v^Bj, 
Jo 



where iV 9 is a Poisson measure with intensity measure A q dsdu, N r is a Poisson 
measure with intensity measure A r dsduQ s {dv Sf ) and -B is a Brownian motion. There 
is existence and uniqueness for this S.D.E. because the Poisson measures N q and N r 
are finite, and because the drift and diffusion coefficients are Lipshitz-continuous 
(see Ikeda-Watanabe dJ). Using Ito's formula and taking expectations, we get, 
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for any ip E C b 2 (K), 

E[<p(V t )] = E[<p(y )] 
rt 



E 



[<p(V a + H«{s, V s ,u)) - <p(V s )]A q du 



ds 



[ip(V s + H r (s, V s ,v ti u))- tp(V s )]A r Q s (dv*)du 



ds 



E[<p'(y a )V a ]d8 + a / E[<p"(V s )]ds 



o 

E[ip(V )} 

rt 

E 
t 



E 



[<p(V, + h)-<p(V.)]q{8,V„dh) 



ds 



[<p(V s +h)- <p(V s )]r(s, V s ,v,,dh)Q s (dv x ) 



ds 



EW{V a )V a ]ds + a E[<p"(V s )]ds 



Setting gf = £(V t ), we thus realize that (gf)t>o solves (A. 4). 

,oo),V 2 { 



If (Qt)t>o and (Rt)t>o are in L^ C ([Q, oc),V2(M.)), then by the same kind of 
arguments as in the uniqueness proof, we have for any t G [0,T], denoting \i t = 



,Q 



9t - 9t 



R 



\h\tv<C / {\hs\tv + \Qs - R s \tv 
Jo 



)ds, 



whence by Gronwall's Lemma, 



sup \ih\tv < Ct \ \Q s -R s Wvds. 
[0,T] Jo 



We consider f e P 2 (M). For t > 0, we set / t ° = f and /, 



fc+i 



rk 

ft . Then we have 



sup \f t 

[0,T] 



fc+i 



ft\ 



TV 



• F \fi - fi-^Tvds. 

Jo 



We clas sicall y conclu de tha t (ft) t>o co nverg es as k tends to infinity to some (ft)t>o 
solving ( |A.3[ ). Using ( |A.l[ ), ( A.2[ ) and ( A.3 | with <p(«) = v 2 , we see that (/ t )t>o € 
Lg= c ([0,c^P 2 (R)). 
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